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Abstract 

Vacuum gravitational fields invariant for a non Abelian Lie al- 

. gebra generated by two Killing fields whose commutator is light-like 

CN| \ are analyzed. It is shown that they represent nonlinear gravitational 

waves obeying to two nonlinear superposition laws. The energy and 

the polarization of this family of waves are explicitely evaluated. 
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In the last years, together with the experimental efforts devoted to the de- 
tection of gravitational waves, there is a strong theoretical activity to describe 
and predict the emission of gravitational waves from astrophysical systems in 
strong field conditions. However, all the experimental devices, laser interfer- 
ometers (Geo-600, Ligo, Virgo, Tama-300) and resonant antennas (Allegro, 
Auriga, Explorer, Nautilus, Niobe), are constructed according to results ob- 
tained by the linearized Einstein field equations, in close analogy with what is 
normally done in Maxwell theory of electromagnetic fields. Indeed, since the 
sources are at enormous distance from the Earth, the amplitude of the waves 
when they reach the detector is so small that it has always been assumed 
that when treating the waves in the Earth's neighborhood the linearized the- 
ory suffices. This assumption is misleading. The nonlinearity of Einstein's 
equations manifests itself in a permanent displacement of the test masses of a 
laser interferometer detector after the passage of a wave train which cannot 
be neglected. Indeed, it has been shown M that every gravitational-wave 
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burst has a nonlinear memory due to the cumulative contribution of the ef- 
fective stress of the gravitational waves themselves and that, in the case of a 
binary coalescence, the nonlinear memory is of the same order of magnitude 
as the maximal amplitude of the dynamical part of the burst. 

On the other hand, the non linearity of the gravitational field is one of its 
most characteristic properties, and it is likely that at least some of the crucial 
properties of the field show themselves only through the nonlinear terms. 
Thus, a great deal of interest is still devoted to exact solutions which more 
easily enable to discriminate between physical and pathological features. 

Starting from the seventy's new powerful mathematical methods have 
been developed to deal with nonlinear evolution equations and their exact 
solutions. Indeed, it has been shown [0 that the vacuum Einstein field 
equations for a metric of the form 

g — f (z, t) (dt 2 — dz 2 ) + tin (z, t) dx 2 + (z, t) dy 2 + 2h\2 {z, t) dxdy 

can be integrated by using a suitable generalization of the Inverse Scattering 
Transform @], leading to solitary wave solutions. 

A geometric inspection of the metric above shows that it is invariant 
under translations along the x,j/-axes, i.e. it admits two Killing fields d x 
and d y closing a Abelian 2 -dimensional Lie algebraQ ^4 2 - Moreover, the 
distribution T>, generated by d x and d y , is 2 -dimensional and the distribution 
T)- 1 , orthogonal to V, is integrable and transversal to T>. 

Thus, it has been natural to consider || the general problem of charac- 
terizing all gravitational fields g invariant for a Lie algebra Q of Killing fields 
such that: 

I. the distribution T> generated by the vector fields of Q is 2-dimensional. 

II. the distribution T> orthogonal to T> is integrable and transversal to 

V. 

The aim of this letter is to study the physical properties of a subfamily of 
such Einstein metrics, better specified below, which has the important char- 
acteristic of representing propagative gravitational fields. The wave character 
of these solutions will be established on a rigorous basis and the problem of 
defining physical observables like the energy flux and the polarization will be 
addressed. 

* The study of metrics invariant for a Abelian 2-dimensional Killing Lie algebra goes 
back to Einstein, Rosen j|, [I| and Kompaneyets recent analysis can be found in |^| [jj. 
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In the following a sub-algebra of the Lie algebra of all Killing fields of 
a metric g will be called Killing algebra of g. Moreover, an integral (2- 
dimensional) submanifold of T> will be called a Killing leaf. 



Invariant metrics. If g is a metric on the space-time and Q2 one of its 
Killing algebras whose generators X, Y satisfy 

[X,Y] = Y, 

then, under the conditions I and II, a coordinate system (x M ), /i = 1,2,3,4, 
adapted to the Killing fields, exists such that 

X = d 3 , Y = e x3 dt, 

and the most general (^-invariant metric has the form 

g = g ij dx l dx j + {\ (a; 4 ) 2 - 2 fix 4 + z/j (dx 3 ) 2 + 
2 (fi - Ax 4 ) dx 3 dx 4 + A (dx 4 ) 2 , z, j = 1, 2 
where g^, A, /x, and f are arbitrary functions of (x^x 2 ). 



Einstein metrics. If the Killing field Y is of light type, i.e. g(Y,Y) = 
0, then the general solution of vacuum Einstein equations, in the adapted 
coordinates (x 1 , x 2 , x 3 , x 4 ), is given || by 

g = 2f((dx 1 ) 2 + (dx 2 ) 2 ) +fi {w(x\x 2 ) -2x 4 )(dx 3 ) 2 + 2dx 3 dx 4 , (1) 

where \x = D$ + B; D, B G TZ, $ is a non constant harmonic function of 
(x 1 ,^ 2 ), / = (V$) 2 ■s/l/j]/ [A, and u> is a solution of the fi-deformed Laplace 
equation: 

Aw + (d 1 In |/i|) dtw + (d 2 In \fi\) «9 2 w = 0. (2) 

Solutions of the /z-deformed Laplace equation will be called fi- harmonic func- 
tions. 

Despite the non-linear nature of general relativity, it is possible to ex- 
hibit two superposition laws for the gravitational fields ([!]). Indeed, with two 
harmonic functions $1 and $2 we can associate three gravitational fields (in 
facts a whole 2-parameters family), that is, g^, g$ 2 and g a <i>i+f>$2i the last 
one, which is associated with a linear combination of $1 and $2, is to be 
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intended as the superposition of the two associated solutions g<s >1 and g^ 2 . 
The second superposition law follows from the linearity of the /z-deformed 
Laplace equation, so that with two /^-harmonic functions Wi and W2 we can 
associate the gravitational fields g Wl , g W2 and g aWl +bw 2 - 

The coordinates (x 3 , x 4 ) on the Killing leaves have a clear geometric mean- 
ing but are of difficult physical interpretation. Fortunately, being the Killing 
leaves flat manifolds |§ it is possible to adopt coordinates diagonalizing the 
metric restricted to the Killing leaves. 

This allows to define harmonic coordinates (x, y, z, t) such that the metrics 
(JJ) take a simpler and physically more interesting form 15). If, for the sake 
of simplicity, we choose /i = 1, then 

x = x 1 

2 

y = x 

z = I [(2x 4 — w (x x , x 2 )) exp (—x 3 ) + exp (x 3 )] 
t = \ [(2x 4 — w (x 1 , x 2 )) exp (—x 3 ) — exp (x 3 )] . 

In these coordinates, the Einstein metrics ([1]) (with /i = 1) read 

g = 2 (V$) 2 [dx 2 + dy 2 ] + dz 2 - dt 2 + dw d\n\z - t\ , (3) 

where t plays the role of time. 

From the physical point of view it is already clear that the Einstein metrics 
(0) represent gravitational waves, the wave-like components, in the almost 
Minkowskian coordinates (x,y,z,t), slowly decreasing when \z — t\ — > oo. 

The choice of harmonic coordinates makes it clear that: 

• when w is constant the Einstein metrics (|^) admit the time-like Killing 
vector field d t , which is orthogonal to space-like hypersurfaces defined 
by t = const, and then represent static gravitational fields. Under the 
further assumption $ = x/v2, they reduce to the Minkowski metric. 

• when w is not constant the Einstein metrics ([3]) represent a disturb 
moving along the z direction on the Killing leaves at light velocity. 

In the following we will assume that w is not constant. The choice fi = 1 
allows to simplify the computational algebra without spoiling the most inter- 
esting features of the metrics like the superposition law and the propagative 
character. Moreover, the form (^) will result particularly manageable for a 
discussion about the definition of spin and energy. Explicit results will be 
illustrated mainly for the case / = const. 
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Wave character of the gravitational field. The second step to attempt 
a physical interpretation of the solutions we are considering is the study of 
their wave character. 

The gravitational fields (Q) satisfy the Zakharov generalization of the 
Zel'manov criterion JT(| which states that a vacuum solution of the Einstein 



equations is a gravitational wave if the corresponding Riemann tensor field R 
is not covariantly constant and its components R^uXa satisfy the hyperbolic 
equation 

g a V ' a V pRfMuXa = KfMuXcn 

where Vp denotes the Levi-Civita covariant derivative and K is a tensor 
field depending at most on first derivatives of the Riemann tensor field. The 
criterion is certainly satisfied if 

g^dadpR^ = 0, (4) 

where d a are now ordinary partial derivatives. 



Concerning the physical meaning of this criterion, it can be shown flO 
that the characteristic hypersurface of the system of equations (f|) is identical 
with the characteristic hypersurface of the Einstein and Maxwell equations 
in curved space-time. Consequently Eqs.(f|) describe the propagation of the 
discontinuities of the second derivatives of the Riemann tensor. This links 
the Zel'manov criterion to the intuitive concept of local wave of curvature. 

In the simple case when the function / is constant the criterion is easily 
verified. In facts, the only independent non vanishing components of the 
Riemann tensor field are 

p d x 2 w d x d y w d y 2 w 

R txzx = 2(2 _ t)2 , R ***u= 2(2 -t)2' Rtyzy = 2(z-ty ] ( ) 

since w is a harmonic function in the plane (x,y), they automatically satisfy 
the condition (|]), the derivatives of w being harmonic functions as well. With 
a little computational effort [[J the result generalizes to /, /x generic. Thus, 
the wave character of the gravitational fields (ff) is established according to 
a rigorous and intrinsic criterion. 

The energy- momentum pseudo-tensor. For / = 1/2, the exact solu- 
tion (H) reduces to 

g = dx 2 + dy 2 + dz 2 — dt 2 + dwd In \z — t\ (6) 
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and has the physically interesting form of a perturbed Minkowski metric, 
g = 7] + h, with rj the Minkowski metric and 

h = dwd\n\z — t\ = [(d x w)dxdz 

+ (dyw)dydz — (d x w)dxdt — (d y w)dydt] (7) 

It is easily verified that the metric (|6|), besides being a exact solution of the 
vacuum Einstein equations, satisfies the linearized equations Dh^ = as 
well. 

To study the energy content and the polarization of gravitational waves 
it is convenient to use the Landau- Lifshitz energy-momentum pseudo-tensor 



H| which has been seen to yield the correct definition of energy for relevant 



cases 



12]. In facts, the energy flux radiated at infinity for an asymptotically 



flat space-time, evaluated with the Landau-Lifshitz pseudo-tensor, has been 



seen to agree with the Bondi flux [13] that is with the energy flux evaluated 



in the exact theory. For gravitational fields (y) it reduces to 
t p« = C{-T» Xa Tl u {f x g^-g^g x ») 

- g^grrZJ^ ~ 9 K V U K^ + g^sTTV^ } (8) 

where the harmonicity condition T^g^ 11 = and the property ff^h^ = 0, 
which in turn implies Y° a = 0, have been used. 

Relevant are the components p^ = Tq of the mixed pseudo-tensor, p^ 
being the 4-momentum density, 

p° = r 2 {t- zy^didjw) 2 + c 2 (d x d yW f] 

+ f- 2 (t-zy 4 C 3 V ■ (\Vw\ 2 Vw) 
p 1 = p 2 = o 

where Cj are positive numerical constants and the harmonicity condition for 
w has been used. 

Concerning the definition of the polarization, the above form for Tq is par- 
ticularly appealing because, apart from a physically irrelevant total deriva- 
tive, the component Tq, representing the energy density, is expressed as a 
sum of square amplitudes. Moreover, the momentum p % = Tq is only non 
vanishing in the direction of the propagation of the wave and it is propor- 
tional to the energy with proportionality constant c = 1 which is the light 
velocity in natural units. 
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The study of the polarization of this wave leads to a non standard result, 
that is h tx = (z — t) 1 d x w and h ty = (z — t)" 1 d y w represent a spin one field. 
This is a consequence of the fact that we are not in the transverse-traceless 
gauge or, roughly speaking, h has only one index in the plane transverse to 
the propagation direction. 

The polarization. The definition and the meaning of spin or polarization 
for a theory, such as general relativity, which is non-linear and possesses a 
much bigger invariance than just the Poincare one, deserve a careful analysis. 

It is well known that the concept of particle together with its degrees of 
freedom like the spin may be only introduced for linear theories (for example 
for the Yang-Mills theories, which are non linear, we need to perform a 
perturbative expansion around the linearized theory). In these theories, when 
Poincare invariant, the particles are classified in terms of the eigenvalues of 
the two Casimir operators of the Poincare group, P 2 and W 2 where P M are 
the translation generators and = ^e flupa P u M pa is the Pauli-Ljubanski 
polarization vector with M pu Lorentz generators. Then, the total angular 
momentum J = L + S is defined in terms of the generators M^ u as J 1 = 
l e °v k Mjf;. The generators P^ and M^ v span the Poincare algebra, XSO(3, 1). 

When P 2 = 0, W 2 — 0, W and P are linearly dependent 

W„ = \P^ 

the constant of proportionality is the helicity. The time component of W is 
W° = ~P ■ 7, so that 

which is the definition of helicity for massless particles like the photons. 

Let us turn now to the gravitational fields represented by Eq. (H). As it 
has been shown, they represent gravitational waves moving at the velocity 
of light, that is, in the would be quantized theory, particles with zero rest 
mass. Thus, if a classification in terms of Poincare group invariants could 
be performed, these waves would belong to the class of unitary (infinite- 
dimensional) representations of the Poincare group characterized by P 2 = 0, 
W 2 = 0. In order for such a classification to be meaningful P 2 and W 2 have 
to be invariants of the theory. This is not the case for general relativity, 
unless we restrict to a subset of transformations selected for example by 
some physical criterion or by experimental constraints. For the solutions of 
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the linearized vacuum Einstein equations the Lorentz transformations are 
selected out by the request that the energy-momentum pseudo-tensor be 
truely a tensor and that the direction of propagation of the waves be preserved 
0- 

Let us see what happens in our case. There exist several equivalent proce- 
dures to evaluate the polarization. One of these consists in looking at the Tq 
component of the Landau-Lifshitz pseudo-tensor, and to see how the metric 
components that appear in Tg transform under an infinitesimal rotation TZ in 
the plane (x, y) transverse to the propagation direction^. This transformation 
preserves the harmonicity condition. The physical components of the metric 
are h tx and h ty and under the infinitesimal rotation TZ in the plane (x, y) 
they transform as components of a vector. Applied to any vector (v 1 ^ 2 ), 
the infinitesimal rotation TZ acts as 

Kv 1 = v 2 , TZv 2 = -v 1 . 

Then, 

TZ 2 v i = -v i z = l,2, 

so that iTZ has eigenvalues ±1. Thus, the components of h^ v which contribute 
to the energy density correspond to spin 1 fields. The reason for this is that 
the propagating part of the metric, in the harmonic coordinates, has a off- 
diagonal form or, put in another way, h^ u has one index in the propagation 
direction and one in the orthogonal plane. As in the linearized theory, of 
the whole diffeomorphisms group just the Lorentz transformations preserve 
the propagation direction and ensure the tensorial character of the energy- 
momentum pseudo-tensor. Moreover, they preserve the harmonic gauge. 

The observable effects of the gravitational waves follow from the study 
of the relative motion of test particles described by the geodesic deviation 
equation. Since the only nonvanishing components of the Riemann tensor 
field R^ au are listed in fl5|), for small velocity and in the weak field approxi- 
mation the geodesic deviation equations reduce to 

d ? V% = 7 -7^ 9 A^ for i = 1, 2, 3 

(z — ty 

TWith respect to the Minkowskian background metric this plane is orthogonal to the 
propagation direction. With respect to the full metric it is transversal to the propagation 
direction and orthogonal only in the limit \z — t\ i— > oo. 
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the vector field V = (0, V 1 , V 2 , V 3 ) representing the space-like separation of 
close geodesies. Thus, it is clear that the deviation depends explicitly on the 
choice of the harmonic function w. A complete analysis will be performed in 
a forthcoming paper by introducing a source for the gravitational waves here 
described. 
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